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Abstract. A theory of electron states for graphene nanoribbons with a smoothly varying width is developed. 
It is demonstrated that the standard adiabatic approximation allowing to neglect the mixing of different 
standing waves is more restrictive for the massless Dirac fermions in graphene than for the conventional 
electron gas. For the case of zigzag boundary conditions, one can expect a well-pronounced conductance 
quantization only for highly excited states. This difference is related to the relativistic Zitterbewegung 
effect in graphene. 
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The experimental discovery of a truly two-dimensional 
allotrope of carbon, graphene p][2] and of the massless 
Dirac character of its electronic energy spectrum [3113] has 
initiated an enormously growing interest in this field (for 
review, see Refs. 15,6 ). One of the most exciting aspects 
of the problem is the hope to develop novel carbon-based 
electronics. Very recently, the experimental realization of 
quantum dots [5] and carbon nanoribbons [3[8] has been 
announced, the former demonstrating single-electron tran- 
sistor (SET) effect [5]. 

The conductance quantization in the ballistic regime [HI 
fl ^ fTT | fT2 l [T3] is one of the most important physical phe- 
nomena determining the functioning of such nanodevices. 
It was considered recently for the case of ideal graphene 
stripe [14] and for the case of confinement due to a smooth 
external electrostatic potential [T5] . The experimental sit- 
uation 5,7 ,8\ corresponds rather to the case of electron 
confinement due to a curvilinear shape of the graphene 
samples than to an external field. The description of the 
penetration of electron waves through constrictions in the 
nanoribbons requires a different theoretical approach. Nu- 
merical calculations of electronic transport in graphene 
nanoribbons demonstrating a very interesting "valley fil- 
ter" effect have been recently published [16]. However, a 
general theoretical analysis of the situation is still absent. 
Here we present an analytical theory of conductance quan- 
tization in graphene nanoribbons based on the adiabatic 
approximation [10, 12. The latter means a separation of 
the electron motion in the directions perpendicular and 
along the stripe. For nonrelativistic electrons the adia- 
batic approximation requires only the smoothness of the 
shape of the stripe boundary and results in the quantiza- 
tion of the conductance. For graphene, the situation turns 
out to be essentially dependent on the boundary condi- 



tions. It will be shown that for the zigzag boundaries this 
theory is essentially different from that for nonrelativistic 
electrons and, in general, there is no reason to expect an 
adiabatic regime and well-pronounced conductance jumps 
for the lowest states of the ribbon. 

The two-component wave function (u, v) for charge 
carriers in graphene with wave vectors close to the K point 
is described by the Dirac equation 



8.8, 
•5- + i-z- I u = ikv, 
ox ay , 



8 
dx 



8y, 



(1) 



where k = E/hvp, E is the electron energy and vf — 10 6 
m/s is the Fermi velocity [3j[4] ; for the other valley K' the 
signs before d/dy are opposite. Let us consider first the 
case of a uniform graphene strip of width L along the y- 
axis, \y\ < L/2. To specify the problem one has to choose 
boundary conditions at the edges p"7] . 

We start with the case of zigzag edges where u(y = —L/2) 
— 0,v(y = L/2) = 0. The energy spectrum is discrete, 
Ej = Kvphj where 



vrj 
L ' 



(2) 



and the wave functions have the form 

u 3 (y) = 77= cos kj (y - L/2) , 

vj (y) = --^= sinfc 3 (y -L/2). 

Consider now the case of a smoothly varying strip 
width, L — > L (x) , \dL/dx\ -C 1. Following a general scheme [101 



(3) 



2 
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1 1 2 j one can try a solution of equation ([T]) as an expansion 

u (x, y) = E c i ( x ) u t ( y ) ' 

3 

v (x, y) = c i ( x ) v f ] (y) ( 4 ) 



where u^ x \v^ xS> are the functions (J3j) with the replacement 
L — > L(x) . The functions ([!]) satisfy by construction the 
boundary conditions. By substituting the expansion ([4]) 
into Eq.JT]), multiplying the first equation by (vj\ and the 
second one by (v,j\ one finds: 



E 

j' l 



dCf_ 

dx 



E 



dcy_ 

dx 



| dVjr 

dx 



| duj> 
dx 



i' 

c f (vj\u f ) , 



j' 



(k 



Cj'(uj\vj>). (5) 



This equation is formally exact. As a first step to the 
adiabatic approximation, one should neglect the terms 



with 



i dv- 



and 



"JtHaT/ \^j\~1x~^ which is justified by the 
smallness of dL/dx, as in the case of nonrelativistic elec- 
trons PIECE]. 

To proceed further we need to calculate the overlap 

L/2 

integrals (<fii\(f)2) = J dy<f>\4>2 for different basis func- 

-L/2 

tions: 



(Vj\Vj>) 
(Uj\Vj>) 



(%' + d ]--j') 
(%' - h-i') 



{Vj>\U 3 ) 



-rferr, j'-j = 2n, W 

"U +3) 



2n + l, 



where n is integer. Substituting Eq.([5]) into Eq.([5]) and ne- 
glecting the nonadiabatic terms with the matrix elements 
of the operator d/dx, we obtain after simple transforma- 
tions: 



dej (x) 
dx 



2i w [k 

IT 



kj> (x)} 



j + f 



(x) 



(7) 



where the sum is over all j' such that j' — j is even. 

Until now we did transformations and approximations 
which are identical to those used in the case of nonrel- 
ativistic electrons. However, we still have a coupling be- 
tween different standing waves so we still cannot prove 
that the electron transmission through the constriction 
is adiabatic. To prove the latter we need one more step, 
namely, a transition from the discrete variable j to real 
one and a replacement of the sums by integrals in the 
right-hand-side of Eq.(j7]): J2'j — — * I $)■■■ where V is 
the symbol of principle value. This step is justified by as- 
suming that kL ^> 1, i.e., only for highly excited states. 



For the low-lying electron standing waves it is difficult to 
see any way to simplify essentially the set of equations ([7]) 
for the coupled states. 

For any function / (z) analytical in the upper (lower) 
complex half-plane one has 



TO 

/ 



dxf (x) 



1 



x\ ± iO 



(8) 



or, equivalently, 



dxf (x) 



x — xi 



±iirf (xi) 



(9) 



Assuming that Cj (x) is analytical in the lower halfplanc 
as a function of complex variable j one obtains, instead of 
Eq.© 

^M = [ k + k j {x)]c^{x). (10) 

Similar, taking into account that c-j (x) is analytical in 
the upper halfplane as a function of complex variable j we 
have 



dc-3 (x) 
dx 



[kj (x) - k] Cj (x) . 



(11) 



At last, differentiating Eq. (fT0|) with respect to x, neglect- 
ing the derivatives of kj (x) due to the smallness of dL/dx 
and taking into account Eq. lfTTj) we find 



d 2 Cj (x) 
dx 2 



+ [k 2 - k) (x)] Cj (x) = 0. 



(12) 



Further analysis completely follows that for the non- 
relativistic case [TO] where k 2 and k 2 (x) play the roles of 
energy and potential, respectively. The potential is qua- 
siclassical for the case of smoothly varying L(x). There- 
fore, the transmission coefficient is very close to one if 
the electron energy exceeds the energy of the jth level 
in the narrowest place of the constriction, and exponen- 
tially small, otherwise. Standard arguments based on the 
Landauer formula [ MT0lfTTiri2l[T3] prove the conductance 
quantization in this situation. 

At the same time, for the lowest energy levels the re- 
placement of sums by integrals in Eq. (fT0)) cannot be jus- 
tified and thus the states with different j's are in general 
coupled even for a smooth constriction ( | dL/dx \ <C 1). 
Therefore electron motion along the stripe is strongly cou- 
pled with that in the perpendicular direction and different 
electron standing waves are essentially entangled. In this 
situation there is no general reason to expect sharp jumps 
and well-defined plateaus in the energy dependence of the 
conductance. This means that the criterion of adiabatic 
approximation is more restrictive for the case of Dirac 
electrons than for the nonrelativistic ones. The formal rea- 
son is an overlap between components of the wave func- 
tions with different pseudospins or, equivalently, between 
hole component of the state j with the electron compo- 
nent of the state j' =/= j. This coupling is a reminiscence 
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of the Zitterbewegung of Dirac electrons determining the 
finite conductivity and anomalous shot noise in graphene 
in the limit of small charge carrier concentration [18U19] . 
Effectively, it works as a kind of intrinsic "disorder" and 
therefore it is not surprising that it destroys the ballis- 
tic regime near the Dirac point. Interestingly, the kinetic 
equation that takes into account the Zitterbewegung ef- 
fects also contains some "P-terms" which are absent in 
the standard Boltzmann equation; these terms become ir- 
relevant for large enough Fermi energy |20j . 

Consider now the case of armchair edges. The bound- 
ary conditions are coupled the components of Dirac spinors 
at K valley u, v with those at K' valley u,v : 

u(-L/2) = u(-L/2), 
v(-L/2) = v(-L/2), 



u(L/2) 
v(L/2) 



e 2 ™u(L/2), 
e 2 ™v(L/2), 



(13) 



where v = 0, ±2/3, depending on the number of rows in 
the strip [T7]. The eigenmodes in that case are just plane 
waves [T7] 



u i (y) = -ivj (y) = 
Uj (y) = -ivj (y) = 



l 



'2L 

1 



exp (ikjy) , 
cxp {—ikjy) 



(14) 

kj = (j + i/Jtt/L, j = 0,±1,±2,... (15) 

A general solution satisfying the boundary conditions (|13|1 
can be probed as 



u(x,y) 


— J]] c j ( x ) ex P ( x ) 

3 


v], 


v(x,y) 


= i bj (x) exp [ikj (x 
j 


)»], 


u(x,y) 


= c i ( x ) exp ^ ik 3 ( x ) 
i 


(L-y)}, 


v(x,y) 


= i bj (x) exp [ikj (x 


){L-y)]. (16) 



Substituting this into the Dirac equation ([I]) one obtains 
exp [ikj (x) v}\-^ + ( kb j - k i c i) + i ~^y°i } = °' 

j 

E, f dbj dkj | 

exp [ikj (x) y] j -± + {kjbj - kcj) + i-^ybj > = 0. 

(17) 

Let us neglect first the nonadiabatic terms propor- 

dk 

tional to in these equations. They will be satisfied 
for sure if all coefficients at the plane waves vanish, which 
is equivalent to the set of equations 



d (Cj + bj) 
dx 

d(bj ~ Cj) 
dx 



+ (k + kj)(bj-Cj) = 0, 

+ (k j -k)(c j + b j ) = 0. (18) 



Differentiating them with respect to x and neglecting, 
again, the derivatives of kj we find the effective Schrodinger 
equation (|12[) and the same equation for bj. Thus, in con- 
trast with the case of zigzag edges, for the armchair edges 
a standard picture of conductance quantization should be 
valid for all states, similar to nonrelativistic electron gas. 

However, there is another problem which makes the 
adiabatic approximation for this case problematic. The 
wave numbers (|15p can depend on x not only due to the 
stripe length but also due to different number of rows in 
the stripe which makes dv/dx a source of sharp random 
potential. It is very difficult to investigate this effect an- 
alytically in the framework of the approach under con- 
sideration. It was argued recently based on numerical re- 
sults and qualitative considerations that this kind of ran- 
domness should be of crucial importance for the graphene 
nanoribbons with the armchair edges [21]. 

It would be very interesting to check experimentally 
the possible difference in the conductance behavior for the 
nanoribbons with zigzag and armchair edges. For the for- 
mer case, the theory predicts essential difference of behav- 
ior at the crossing of low-lying and highly excited energy 
levels in the quantum point contact situation, that is, for 
a narrow constriction of the graphene nanoribbons. 
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